Singlet structure function gi at small x and small 
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Explicit expressions for the singlet gi at small x and small are obtained with the total resum- 
mation of the leading logarithmic contributions. It is shown that gi practically does not depend on 
X in this kinematic region. In contrast, it would be interesting to investigate its dependence on the 
invariant energy 2pq because, being g\ positive at small 2pq, it can turn negative at greater values 
of this variable. The position of the turning point is sensitive to the ratio between the initial quark 
and gluon densities, so its experimental detection would enable to estimate this ratio. 
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I. INTRODUCTION 

The Standard Approach (SA) for the theoretical description of the spin structure function gi at large and small 
X is based on the Altarelli-Parisi or DGLAP Q^— evolution equations Q complemented with global fits for the 
initial parton densities. Originally, SA was suggested for describing the region of large x but later it has been applied 
for investigating the polarized DIS at small x as well. As SA neglects the total resummation of leading ln(l/a:), 
which becomes necessary at small x, the singular (~ factors are introduced in the fits for the initial parton 

densities. As it was shown in Refs. 0, 0, such factors act as the leading singularities^ in the Mellin space. They 
ensure the steep rise of 51 at a; <C 1 and indeed mimic the impact of the total resummation of leading ln(l/x) terms. 
Alternatively, when the total resummation is taken into account, those singular factors become unnecessary, so the 
initial parton densities can be fitted with much simpler expressions. The total resummation of ln(l/a;) contributions 
to the anomalous dimensions and the coefficient functions of the singlet component of gi was done in Ref. 5] in the 
Double- Logarithmic Approximation under the assumption of as fixed at an unknown scale. More precise results 
including the running as effects were obtained in Ref. ifij. 

In the present paper, we extend the results of Ref. 0] to consider the small- x behavior of the singlet gi in more 
detail. In particular, we give a special attention to the kinematic region where not only x but also are small. 
On one hand, this kinematics has been investigated experimentally by the COMPASS collaboration, see Ref. 0. On 
the other hand, the region of small is clearly beyond the reach of SA. We show that in this kinematics gi can be 
practically independent of x even for x We obtain that 51, being positive at small values of the invariant energy 
2pq, can turn negative when 2pq increases. The position of the turning point is sensitive to the ratio between the 
initial quark and gluon densities. Then we also show that, in spite of the presence of large factors providing gi with 
the Regge behavior at small a;, the interplay between initial quark and gluon densities might keep gi at small x close 
to zero even at small x, regardless of values of Q^. 

The paper is organized as follows: in Sect. 2 we remind and explain the basic formulae for gi singlet obtained in 
Ref. These formulae include the total resummaion of the leading logarithms of x. In our approach, the coefficient 
functions for gi are expressed through new anomalous dimensions. Explicit expressions for them are presented in 
Sect. 3. We focus on gi at small in Sect. 4. As our approach is perturbative, we are interested in minimizing the 
influence of non-perturbative contributions. To his aim we introduce an optimal mass scale in Sec. 5. The asymptotics 
of gi at small x is considered in Sect. 6. Suggestions for new simple fits for the initial parton densities at small are 
briefly discussed in Sec. 7. Sect. 8 contains our numerical results, and flnally Sect. 9 is for the concluding remarks. 



They are simple poles whereas the total resummation leads to the leading singularity as the square root branch point. 
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II. EXPRESSIONS FOR gi AT LARGE 

The singlet structure function gi at small x was studied in Ref. |6||. According to it, gi can be represented in the 
form of the Mellin integral: 



where < > stands for the sum of electric charges: < >= 10/9 for n/ = 4, i/ = ln(Q //i ), with /i being the 
starting point of the Q^- evolution, 8q is the initial averaged quark density: < e^^> 8q — e\bu + e\bd + ... whereas 3g 
is the initial giuon density. 

The other ingredients of the integrand in Eq. are expressed in terms of the anomalous dimensions Hik , with 
i,k = q,g. The exponents f^(±) and coefficient functions Cq^^ are: 

n^±)^l[Hqq+Hqq±R] , (2) 
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with -q = ln(/i2/AQ(^^) and b = (33 — 2n^)/(127r). The additional factor ^ — jmJ^^ ™ coefficients Cg^'' is the 

small-w estimate for the quark box which relates the initial gluons to the electromagnetic current. A'{lu) stands for 
the QCD coupling in the box in the Mellin space. 



III. ANOMALOUS DIMENSIONS 

The anomalous dimensions Hik obey the following system of equations: 
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with the Born contributions Uik defined as follows 
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A' is given by the Eq. (O and 
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-if+TT"^ Jo (p + 77)2+7r2 

is the Melhn representation of the QCD running couphng as involved in the quark-gluon ladder, with the proper 
account of its analytic properties. In Eq.® we use the standard notations for Cp = {N^ — 1)/{2N) = 4/3 and iV = 3. 
Finally, 
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is the factor that accounts for non-ladder diagrams. 
The solution to Eq. © is 
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where 



Z = - 2(6,, + h,g)) + ^/(^2 _ 2(6,, + 6„))2 - 4(6„ - 6„)2 - 166„6,g 
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IV. EXPRESSIONS FOR THE SINGLET gi AT SMALL 

Eq. states that gi does not depend on when ^ ji^ . In this case gi depends only on z = ji^ /{2pq): 
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Eq. ^ was obtained for > ji^ and cannot be used for studying the -dependence of gi at < ^ . However, it 
would be interesting to extend our approach to this region. In order to do so, we suggest to modify Eq. replacing 
Q2 by (Q2 -(- ^2^^ Although such a shift takes us out of the logarithmic accuracy we have always kept in our approach, 
it looks quite reasonable and natural. Indeed, let us in the first place consider a contribution of a ladder Feynman 
graph at a; ^ 1 with the DL accuracy. The graph includes the quark and gluon rungs. Integrations of the quark 
rungs are infrared-stable, being regulated with the quark mass m,. On the contrary, integrations of the gluon rungs 
are IR-divergent, so they must be regulated. The standard way of the IR-regulating in QED and QCD is providing 
gluons with a mass [i which acts as an IR cut-off. It is also convenient to choose \x ^ m, and replace m, by \x in the 
quark propagators as was first suggested in Ref . @ . After that m, can be dropped. Now both gluon and quark rungs 
of the ladder are IR-stable and [i -dependent. The simplification of the spin structure can be done with the standard 
means (see e.g. the review "9]). It is appropriate to use the standard Sudakov variables for integrations over momenta 
ki of ladder virtual quarks and gluons: ki — ai{p — {m'^/2pq)) + l3i{q + xp) + k±. After that the DL contribution of a 
ladder graph with n rungs is proportional to the integral J„: 
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where the rungs are numbered from the bottom to the top of the ladder. We have used the notation w = 2pq. As we 
consider a; ^ 1, we neglected the term —wxan coming from the representation 2gfc„ = ?ii/3„ — wxan in the argument 
of the first 5 -function and similar terms in the other 5 -functions. Eq. 1)16(1 manifests that the -dependence in J„ 
at a: ^ 1 is given by the term + fi^ in the first S -function only. Neither accounting for non-ladder graphs nor 
accounting for single logarithms, including the running as effects, change this situation. So, the replacement 

is confirmed by the analysis of the structure of the Feynman graphs, though as the replacement is beyond the 
logarithmic accuracy, it can be called model-dependent. 

Further, DGLAP exploits the -evolution, with being the upper limit of integrations over fc^ ^ (r = 1,2, ..), 
so the first loop DGLAP double-logarithmic contribution jDGLAP _ ln(Q^//i^). This contribution is large only when 
3> /i^. In contrast, we use the evolution with respect to /i^, the integrations over fc^ ^ run up to w instead of , 
so our approach is not restricted by the region of large . For example, when n = 1, Eq. (|16() yields 

Ji = ln(w/^2) ^ (^8^ 
so Ji does not depend on at all. The -dependence appears in the next loops. In particular, when n — 2, 

J2 = -(1/2) ln2(w//?) HQ^w) + (1/6) In'iQ^/w) , (19) 

so it depends on through Q^ + jj.^. It agrees with Eq. (|17|) . In contrast to DGLAP, double logarithms in Eqs. (|18I19|) 
do not disappear when ^ 0. 

Using Eq. H17|l makes possible to rewrite Eq. Q in the form convenient equally for large and small Q^: 

-<-«"'^^/:£(i^)i(-^^'<-K^)"'"-j-'M(2!±i^)"'->, ,0, 

-^(-s^'<-')(^)'"'--i-'<-'K^)°'-)H- 

Eq. (pn|) coincides with Eq. (QJ when ^ and also reproduces Eq. ((T3|l when = 0. Eq. shows that 
the X- and Q^- dependence of gi are getting weaker with decreasing so that gi at <C /i^ depends rather on 
z = fi'^/{2pq) than on x or Q^. Eqs. (|1I2U|) describe the leading Q"^ -dependence of gi at :$> Similarly, Eq. 
describes the leading -dependence at <C /i^: although logarithms of {{Q^ + ji^)/ ^"^^ are small here, they are 
multiplied by leading, double logarithms of l/z contrary to other, unaccounted -terms. In what follows we will 
not discuss the Q^- dependence of gi at small in detail. Instead, we focus on investigating gi at ^ where gi 
is given by any of Eqs. H15I20I) . According to Eqs. 11I15I20|I . the total resummation of double-logarithms for gi makes 
it depend on the value and the way of introducing the cut-off /i. This dependence will vanish when the probabilities 
to find a polarized quark and gluon are calculated and used in expressions for gi instead of the phenomenological 
initial densities 5q and 5g. 



V. OPTIMAL SCALE FOR m 

In order to estimate we discuss below the restrictions for it. From 

as fc^ ^ ^QCD > is^, we obtain the obvious restriction for the value of /x: 

» ^QCD ■ (22) 

Then, the DL contributions from ladder quark rungs are infrared- stable, with logarithms there containing masses 
niq of the involved quarks in denominators. In order to calculate ladder fermion graphs with Infra-Red Evolution 
Equations these logs should be regulated with the infrared cut-off fi. It brings the second restriction for /i: 



fj, > m, 



(23) 
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Basically, there are no other restrictions for /i. However, some additional information of /i comes from the small- a; 
asymptotics of 171 . In Ref . &l it was shown that 

91 - (l/xr (24) 

when X — > 0. It turned out that ujq depends on 77 = ln(/i^/AQp^), in such a way that loq is maximal at 77 = 775 ~ 7.5 
with 



wo(?7s) = As« 0.86. (25) 

We have called A5 the intercept of the singlet gi. This value is in agreement with the analysis of experimental 
data^3|- Assuming Aqcd = 0.1 GeV leads to the estimate 

MS = AQCoe^-^'' « 5.5 GeV . (26) 

On the other hand, from physical considerations, the intercept A5 should be a constant and should not depend on 
/X. This dependence is the artefact of our approach: we account for perturbativc contributions to the asymptotics 
and leave out a possible impact of non-perturbative ones. Taken together, the perturbative and non-perturbative 
contributions would make loq to be fi -independent. We suggested in Ref. [a| that non-perturbative contributions to 
ujQ happened to be minimal at 77 = 77s so in order to minimize the impact of (the basically unknown) non-perturbative 
contributions on (71, we should fix /i = /ug = 5.5 GeV. We call the Optimal Scale for the singlet gi. We expect 
that choosing this scale for /i would bring a better agreement between experimental data and our formulae than other 
values of /i. We also suggest that the initial parton densities can be fitted mostly simply when 7/ is fixed at the 
Optimal Scale. It is worth of mentioning that in Refs. [TH the Optimal Scale fj,NS for the non-singiet component of 
gi is 5 times smaller: fj,Ns — 1 GeV. 



VI. SMALL- X ASYMPTOTICS OF gi 

Before performing numerical analysis of Eq. I|2U|) . it could be instructive to consider its small- a; asymptotics. This 
asymptotics is different for small and large and when a; + z — > 0, we obtain 
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where the intercept A5 is given by Eq. H25() are taken at the intercept point uj = As-, A5 w 0.86, 
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with all hij and their ^-derivatives in Eq. (|28|l are taken at the intercept point to — A5. The initial parton densities 
5q{uo) and Sg{uj) are also fixed at = A5. 

When z and ^ , Eq. (|27l) turns (we drop here the unessential overall factor) to 



gi^5(A5)<5<7(As)(i)^'/ln'/'(lA) 



(30) 



with 



S{As) = -1 - 0.064 5g{As)/5q{As) 



(31) 



Eqs. H30I31|I show that the aymptotics of gi does not depend on x and in the small-Q^ region and the sign of gi 
is determined by S{As). There can be three options : 
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A. Large and negative Sg: S{As) > 0. 

When the initial gluon density is negative and large so that 

6g < -IbMSq , 

the asymptotics of gi is positive. It is known that > at large z where it is given by its Born expression 
if Sq and 6g are related by Eq. ()32|l . gi is positive in the whole range of z. 

B. Positive or small and negative Sg: S{As) < 0. 

On the contrary, when 

Sg > -15.64(5, , (33) 
gi, being positive at large x, should pass through the zero value and changes its sign at asymptotically small z. 

C. Fine tuning: 5(As) = 

Finally, there might be a strong correlation between Sq and Sg : 

Sg = IbMSq (34) 

when z — !■ 0. In this case gi is positive at large x and then gi — s- in spite of the large power-like factor (1/z)'^^ in 
Eq. jSOl)- 

VII. FITS FOR INITIAL PARTON DENSITIES 

In the Standard Approach, the initial parton densities Sq{x), Sg{x) are fitted from experimental data at a; ~ 1 and 
= fi^ K 1 GeV^. Then they are evolved with the anomalous dimensions into the region of large : ^ /i^ 
and finally evolved with the coefficient functions into the region of a; ^ 1. As the coefficient functions in the SA do 
not include the total resummation of In a; and therefore cannot provide gi with the steep rise at small x, this role 
is assigned, in the standard fitsQ, to the singular factors which mimic the resummation. In other words, the 
impact of the NLO terms in the DGLAP coefficient functions on the small- a; behavior of gi is actually negligibly 
small compared to the impact of the fit. When the resummation is accounted for, the singular factors can be dropped 
out and the fits can be simplified down to expressions ^ Nq^g{l + Cq^gx). Obviously, the straightforward evolution of 
the fits backwards, to the region of <C ^Ji'^ is beyond SA. We suggest that the analyses of the large Q^- and small 
Q^- experimental data would be more consistent when the argument x in the new fits is replaced by (z -f x). This 
argument « x at large and w z at small . It means that at small the fits should depend on 2pq only. We 
suggest that the fits for Sq(z + x), Sg{z + x) can be chosen as linear forms ~ Nq g{l + aq^g{z + x)), with parameters 
Nq^g, Qq^g fitted from experimental data either at small or large Q^. However, in the present paper we do not plan 
to study new fits in detail. Instead, we assume that after the total resummations of In a; has been accounted for, one 
can approximate the initial parton densities by constants: 

Sq^Nq, Sg^Ng . (35) 

The analysis of the data on the non-singlet gi shows that Nq should be positive whereas the sign of Ng is basically 
unknown. 

VIII. NUMERICAL RESULTS FOR gi AT = 

Both Eq. (|20() and the asymptotic expression Eq. (|30() imply that the x- dependence of gi at small should be 
weak even for very small x. As a matter of fact, gi in this region depends on 2pq only, so a plot of gi vs x should 
reveal a very flat behavior, with the magnitude, (/i(z), depending on the interplay between Sq and Sg at different 



(32) 
Therefore, 
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FIG. 1: Gi evolution with decreasing z = /2[pq) for different values of ratio r — 5g/5q: curve 1 - for r = 0, curve 2 - for 
r — —5 , curve 3 -for r = —8 and curve 4 -for r = — 15. 



z — fi^ /{2pq). Presuming that 6q > Q and approximating 5q and 5g by the constants Nq^g, we rewrite Eq. (|20|) at 

g2 = as 

gi(z) = (<e2>/2)7V,Gi(z) (36) 

and calculate Gi numerically. The results for different values of the ratio r — Ng/Nq, Gi are plotted in Fig. 
When the gluon density is neglected, i.e. Ng = (curve 1), Gi being positive at x '--^ 1, is getting negative very soon, 
at z < 0.5 and falls fast with decreasing z. When Ng/Nq = —5 (curve 2), Gi remains positive and not large until 
z ^ 10^^, turns negative at z '--^ 0.03 and falls afterwards rapidly with decreasing z . This turning point where Gi 
changes its sign is very sensitive to the magnitude of the ratio r . For instance, at Ng/Nq = —8 (curve 3), Gi passes 
through zero at z ~ 10"'^. When Ng/Nq < —10, Gi is positive at any experimentally reachable z (curve 4) . Therefore, 
the experimental measurement of the turning point would allow to draw conclusions on the interplay between the 
initial quark and gluon densities. 

IX. CONCLUSION 

We have shown that the study of 51 at small-Q^ could be as interesting as in the large-Q^ region. Eq. H20II describes 
the singlet gi at small x and arbitrary values of Q^, generalizing both Standard Approach and our previous results. 
It accounts for the total resummation of the leading logarithms of x and z (z = ji^ /2pq). It makes possible to simplify 
the fits for the initial parton densities. In general, gi includes both perturbative and non-perturbative contributions. 
In order to minimize the impact of the latter, we have introduced an Optimal Scale for /z. In the small- region, 
Eq. (|20|l predicts that gi essentially depends on 2pq only and practically does not depend on and x even at a; <C 1, 
making the investigation of the x -dependence uninteresting. On the contrary, the study of the z -dependence of gi at 
small would be useful. Indeed, the sign of gi is positive at z close to 1 and can remain positive or become negative 
at smaller z, depending on the ratio between 5g and 5q. Our numerical results are plotted in Fig. 1. The position of 
this point is sensitive to the ratio Sg/Sq, so the experimental measurement of this point would enable to estimate the 
impact of Sg. 
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^ We remind that Eq. 1201 is becoming unrehable when 2 ~ 1 and should be modified according to the prescription of Ref. in order to 
describe this region. 
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